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Abstract 

We discuss a model nonfrustrated magnetoelectric in which strong enough magnetoelectric cou- 
pling produces incommensurate magnetic order leading to ferroelectricity. Properties of the mag- 
netoelectric in magnetic field directed perpendicular to wave vector describing the spin helix are 
considered in detail. Analysis of classical energy shows that in contrast to naive expectation the 
onset of ferroelectricity takes place at a field H^-i that is lower than the saturation field -ffc2- One 
has Hci = Hc2 at strong enough magnetoelectric coupling. We show that at H = the ferro- 
electricity appears at T = TpE < Tjy. Qualitative discussion of phase diagram \n H — T plane is 
presented within mean field approach. 

PACS numbers: 75.80.+q, 71.70.Ej, 77.80.-e 
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I. INTRODUCTION 



In the past several years, there has been a revival of interest in magnetic ferroelectrics in 
which magnetic and ferroelectric orders coexist (magnetoelectrics) .-i^ Of particular interest 
now are systems in which ferroelectric and spiral magnetic order occur simultaneously due 
to recognition of the role that such materials might play in fabricating novel magnetoelec- 
tric (ME) devices. - A number of such compounds have been obtained recently: i^MnOs 
with R =Gd, Tb, Dy;^^ i^MusOs with R =Ho, Y, Tb, Dy^^ NiaVaOs;-^^^ spinel oxides 
i?Cr204 with R =Co, Fe, Mn;— MnW04,— etc. In the majority of these materials paraelec- 
tric phase with coUinear spin structure (sinusoidal spin density wave) appears below Neel 
temperature T/y. Upon further cooling transition to the phase takes place aX T = TpE < T/v 
in which ferroelectric order coexists with incommensurate elliptical (conical in i?Cr204) mag- 
netic spiral. Such ferroelectric phases are stable down to very small temperature in some of 
them (e.g. i^MnOs and i?Cr204) while in others there are transitions to collinear paraelectric 
phases below TpE- All experiments point to key role of the noncoUinear spin configurations 
induced by frustrated exchange interactions in producing the electric polarization.- Due to 
frustration T/y and TpE are quite small in all compounds found by now with only one excep- 
tion, Bao.sSri.sZnaFeiaOaa,— in which TpE is greater than room temperature. A "giant" ME 
effect is observed in these materials lying in very high sensitivity of the electric polarization 
to the magnetic field: spin-flop transition in magnetic field is accompanied by rotation of 
the polarization by 90° and by anomaly in dielectric constant. Meantime the value of elec- 
tric polarization was found to be two-three orders smaller than in the typical ferroelectrics 
and there is no such great influence of electric field on the magnetic properties indicating 
smallness of ME coupling in these compounds. Then, many efforts are made now to find 
materials with stronger ME coupling and with higher transition temperatures which can be 
used in practice. 

Phenomenological treatment of the mechanism of magnetoelectric coupling has been pro- 
posed basing on Landau expansion and symmetry consideration.-»i^'i^»i^ A microscopic mech- 
anism of the ferroelectricity of magnetic origin has recently been proposed in Ref.-iS, which 
is based on the idea that spin current js oc [Sj x Sj] is induced between the noncollinear 
spins that leads to the electric moment P oc [eij x j^], where ejj is the unit vector connected 
spins i and j. This result can be regarded as an inverse effect of the Dzyaloshinskii-Moriya 
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interaction. As a result one can write the effective ME interaction in the following form:— 

Vme = /5[U X e,,] ■ [S, X S,] , (1) 

where U stands for the corresponding ligand displacement. Taking into account the elastic 
energy, 7f/^/2, one finds that the proposed mechanism can lead to ferroelectricity as soon 
as noncoUinear spin structure exists. 

Meantime ME coupling ([T]) can produce a spiral incommensurate magnetic order and 
electric moment even without frustration if (3 is large enough. Really, let us consider two 
spins and take into account the direct exchange coupling J between them, ME interaction 
([1]) and elastic energy 7?7^/2. Minimization of the total energy with respect to U and cf), 
angle between spins, gives U = {P / ■y) S"^ sin (j) and sin0cos0 = 7J/(/?S')^ sin 0. Meantime 
ME coupling constant in TbMnOs, one of the best magnetoelectric of this type by now, is 
estimatec^^ to be /3 ~ 1 meV/A. Taking into account that the characteristic value of 7 is 
10^ meV/A^ and J ~ 1 meV we have from the above estimations sin cos ~ 10^ sin that 
leads to = and indicates that in magnetoelectric compounds have been found nowadays 
frustration is really indispensable for appearance of noncoUinear magnetic order leading to 
ferroelectricity. As we mention above, every efforts are made now to find incommensurate 
magnetoelectrics with larger ferroelectric moment. But strong ME interaction can produce 
spiral incommensurate magnetic order and ferroelectric moment without frustration. Thus, 
it is seen from the above consideration of two spins that if (3 were about 30 times larger 
than in TbMnOs a nonzero solution for appears. Moreover, one can expect that among 
magnetoelectric compounds with strong enough (3 to be obtained (as we hope) nonfrustrated 
ones would have larger transition temperatures than frustrated ones, other things being 
equal. 

Thus, it would be instructive to discuss nonfrustrated magnets with strong ME coupling of 
the form ([1]). Such magnet is considered recently in Ref.-i^, where collective magnetoelectric 
modes are discussed. In the present paper we discuss a model similar to that of Ref.— 
focusing on its properties in magnetic field directed perpendicular to wave vector describing 
the spin helix. Analysis of classical energy presented in Sec. [Ill shows that in contrast to 
naive expectation the onset of ferroelectricity takes place at a field Hd that is lower than 
the saturation field Hc2- One has Hd = Hc2 at strong enough (3. We show that at H = the 
ferroelectricity appears at T = TpE < T^- Qualitative discussion of phase diagram in H — T 
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FIG. 1: Projections of spins on xy and xz planes are shown in the non-colhnear phase of mag- 
netoelectric discussed. Magnetic field H, electrical polarization P and canting angle of spins in 
magnetic field are depicted. 



plane is presented within mean field approach in Sec. IIIII Sec. IIVI contains our conclusion. 



II. CLASSICAL ENERGY 

We discuss a magnetoelectric with ferromagnetic interaction J![J^ in xy plane and antifer- 
romagnetic one J^j along z axis with ME coupling of the form ([T]) which Hamiltonian has 
the form 

^ = ^ E "^'S^S, + /? J][Uh, X e.] . [Sr. X Sn^+eJ + | E + ^ E 

(2) 

where and j) denote nearest neighbors along z axis and in xy plane, respectively, 
e^, is the unit vector along x axis, lattice constant is taken to be equal to unity, (3 and 7 
are positive constants and the last term is the Zeeman energy in the field directed along z 
axis (see Fig. [1]). At if = spins lie in xy plane and their rotation describes by wave vector 
q = (g, 0, 0). There is a uniform displacement along y axis Uj = U = (0, t/, 0). The electric 
polarization of the sample P is proportional to A^U, where N is the number of spins in the 
lattice. When H ^ the spins cant in opposition to the field direction and makes an angle 
6 < 7t/2 with z axis (see Fig. [T]). Magnetoelectrics with ferromagnetic exchange along z 
axis or antiferromagnetic exchange in xy plane can be considered on the equal footing. We 
discuss corresponding results qualitatively in Sec. HVl 

To find q, U and 9 one has to minimize the classical energy according to U , q and 9 that 
has the form 

^ = 2rS'^ cos^ 9 - r^S^ {cos^ 9 + sin^ 9 cos q) - PS^U sin^ ^ sin g + ^U^ - HS cos 9, (3) 
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where 5" is the spin value. Eq. ([3]) has two solutions: (i) that with collinear spin structure, 
q = U = 0, and (ii) that with spiral spin structure, g 7^ and U 0. The last one has the 
form 

U = sin^ ^ sin g, (4a) 

7 

j^y = cosgsin^6', (4b) 

2/525-3 

H = — — cos e (sin^ e + 2j, - j^y) , (4c) 
where two dimensionless constants are introduced 

Stability conditions of the solutions are determined from the demand of positive definiteness 
of the bilinear form d'^E/{dxdy), where x,y = U,q,6. In particular the stability criteria of 
dl]) are given by 

q 7^ 0,^ (6a) 

cos^ < yiZ^i+^Zi. (6b) 

It is seen from Eq. that this solution exists if 

3.y < 1. (7) 

We assume below that condition ([7]) holds. The second solution of Eq. gives collinear 
spin structure: 

U = q = 0, (8a) 

f H/Hc2, ifH< H^2, 
cos9 = { ' ' - ' (8b) 

[1, iiH>H,2, 

where Hc2 = 4:3 J^. Solution is stable at large enough fields so that the following condition 
satisfies: 



cose>^l-3^y. (9) 

One concludes from Eqs. ([7]) and (jH]) that if ([7]) does not hold only collinear spin structure 
exists. In contrast, when ([7]) is satisfied, the collinear solution is stable only a.t H > H^°, 
where 



H^,° = 4Sr,/T^y (10) 
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FIG. 2: Sketch of the right part of Eq. as a function of cos^. Values of H and 9 are given by 
Eqs. (fT3]) and (fT2]) . respectively. 

that is found using Eq. ^p) and assuming the equaUty in Eq. ([9]). The angle 0^° correspond- 
ing to the field Hl° is given by 



Notice that at = Eq. (Hd) gives g = 0. 

Let us turn to the transition between spiral and collinear configurations. At if = 
and 9 = Ti/2 spiral configuration is realized. As it is clear from Eq. (j3t) the angle 9 is an 
ambiguous function of H. A sketch of the right part of Eq. (jH^) is presented in Fig. [2l It is 
seen that there is a maximum at ^ = where 



If the right part in Eq. (1121) is larger than unity, it indicates that ^ = in consideration 
presented below and the regime with 9 > 9 remains only. According to Eq. ( ]6b[) one should 
take those solutions of Eq. (jH^) which has 9 > 9. The value of magnetic field corresponding 
to the maximum in Fig. [2] is given by 



It can be easily shown using Eqs. ( JTOl) . (|T3l) and Cauchy's inequality that H > H^". Then, 
analysis shows that the type of phase transition is determined by the value of the angle 
9^° at which inequality iQ turns into equality. Two regimes are possible at which we have 
continuous and discontinuous transitions, respectively: 9^° < 9 and 9^° > 9. It can be easily 
shown using Eqs. ( fTTI) and ( fT2|) that 91° < {>)9 is equivalent to jxy + jz < (>)1- Let us 
discuss these two regimes separately. 



COS^^l° = - Jxy 



(11) 




(12) 




(13) 
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A. Continuous transition 



The transition is continuous if 91° > 6, i.e., if 

jxy+Jz>l- (14) 

The angle 6 lowers as the field rises and the spiral solution turns into collinear one at 
H = H^°: q and U reduce gradually to zero as H approaches H'^° and 9 approaches at 
H = Hl° we have 9 = 9^°, q = U = 0, the stability criterion of the spiral solution fl^al) ceases 
to hold and the collinear solution ([8]) becomes stable (criterion begin to hold). All spins 
become parallel the field direction at H = Hc2- As a result we obtain phase diagram shown 
in Fig. [3](a) corresponding to the line T = 0. 



B. Discontinuous transition 



The transition is discontinuous if 9^° < 9, i.e., if 



jxy+jz<^- (15) 

This regime corresponds to larger spin- lattice coupling (larger P) than that discussed above. 
In this case the angle 9ci can not be reached gradually because, in particular, the spiral 
solution is unstable aX 9 < 9 > 9^°. Thus the transition is of the first order in this case. 
When H reaches if^f , q does not turn into zero and the spiral solution remains stable. At 
the same time the collinear solution is also stable aX H > but the energy of spiral 
solution is lower than that of collinear one aX H = H'^°. As one increases the magnetic field 
further, the ground state energies of these two solutions comes together and the transition 
takes place when they become equal. Corresponding field can be greater or lower than Hc2- 
In the first case all spins in the collinear phase are parallel to the field (^f° = 0) whereas 
in the second case 9i° ^ and all spins becomes parallel each other at if > Hc2 only. It is 
easy to find in the second scenario using Eqs. ([3]), (jlj) and ([8]) for the critical field 

Hj:," = 2sr ^~^''i^^\ (16) 



Jz 
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The transition is accompanied by the drop of the angle 9 from 9^^ given by 

cos9i° = ^, (17) 



to 9^,1 



;° given by Eq. (fTTl) and the wave vector q of the spiral switches from 



cos 



Jo _ 



Jxy 



(18) 



to g = 0. Then we lead to the part of phase diagram m H — T plane shown in Fig. [3t^b) 
corresponding to the line T = 0. 
One obtain that iff ° < Hc2 if 



If (3 is as large as this criterion does not hold, there is only one critical field He and we 
lead to the line T = on the phase diagram shown in Fig. El^c). Expression for H^. is quite 
complex and we do not present it here. 

Let us discuss now the phase diagrams at T > 0. 

III. FINITE TEMPERATURES 

Let us find the equation for the line of phase transitions between coUinear and spiral 
phases within the mean-field approximation (MFA). The energy ([3]) is a function of T. 
Working in MFA we imply that the spin value is reduced by thermal fiuctuations: 



where T^v = 2S'^{J'' + 23="^) within MFA, 5 > 1 and T < T^. For simplicity we neglect 
dependence of constants [3, 7 and J on T in Eq. (IHl). As soon as the spin value depends 
on T, all quantities depending on 5* arising in the above discussion are also functions of T. 
In particular, using Eqs. (ITU]) and (1161) one can find the line of phase transitions in H — T 
plane between spiral and collinear phases. Meantime this line is different in the case of 
second and first order phase transitions at T = 0, i.e., at jxyiT = 0) + jz(T = 0) > 1 and 
ixyiT = 0) + jz(T = 0) < 1, that should be discussed separately. 

A. Second order phase transition at T = (i.e., jxy{T = 0) + jz{T = 0) > 1) 

Evidently, the line of phase transition started at [H = 0, T/v) should lead to the point 
{Hc2,T = 0). It is also clear that the line of phase transition between non-collinear and 
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FIG. 3: Phase diagram in the H — T plane of the magnetoelectric discussed in some hmiting cases 
indicated within frame in each plane (constants j are given by Eqs. ^ and j(0) = j{T = 0)). 
There is an incommensurate magnetic order with wave vector q and electric moment P in phase 
I. Spins are parallel each other and canted by a finite angle to the magnetic field in paraelectric 
phase II. All spins are parallel to the magnetic field and P = in phase III. Solid and dashed lines 
denote lines of second and first order phase transitions, respectively. One goes on successively from 
plane (a) to (b) and (c) increasing (3. (a) The line of phase transitions between I and II is given 
by Eq. (I2ip within MFA. (b) Transitions between I and II is of the second and of the first order 
above and below T*, respectively, where T* is given by Eq. (j22p . The line of phase transitions in 
this case is given within MFA by Eqs. (j2ip and (j23p above and below T* , respectively, (c) There is 
only one critical field He at T = in this regime. Temperatures T* and T** are given by Eqs. (j22p 
and (pl|) . respectively. 



coUinear phases should start at {Hci,0) and end at {0,Tfe), where Tfe < T^. Using 
Eqs. f fTOj) and fl20|) one can easily find the equation on H^°{T): 

T [Him? - [Hli{T)f 



N 



2{4SJ 



z\2 



(21) 
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that gives a parabola in H — T plane (see Fig. [3](a)). In particular, we have from Eq. (12T|) 
TpE = Tn[^ ~ jxy{0)]/2. Then the demand T ^ T/v implies that jxy{0) ~ 1. 



B. First order phase transition at T = (i.e., jxy{T = 0) + jz{T = 0) < 1) 

The phase diagram is different depending on whether condition ( IT9|) holds or does not 
hold at T = 0. Let us start with the first case. 



1- V1^ + VjM>i 

It is seen from Eqs. ([5]) that thermal fluctuations increase jxy and jz- Then, above 
a certain temperature T* the sum jxyiT) + j^(T) becomes larger than unity. Hence, at 
T > T* and T < T*, where 

T* = T^^^J^M^lM (22) 

that is found from the condition jxy{T*) + jz{T*) = 1, the transition is of the second and 
first orders, respectively. The line of the second order phase transitions at T > T* is given 
by Eq. (!2T]) . The line of the first order phase transitions at T < T* can be found using 
Eq. flTB]) with the result 

TV = \s3' (2^' 
Then, we lead to plane (b) in Fig. [31 



^- A/Mo) + \/iI(o) < 1 

Phase diagram in this case is presented in Fig. [3]^c). The temperature T** is found from 
the condition ^/ jxyiT**) + ^/JjT**) = 1 with the result 

= Tm [l -slh^)-VJM^ ■ (24) 

At T < T** there are transitions from spiral phase to collinear one with all spins are along the 
field direction (phase III). In contrast, at T > T** there is first a transition to collinear phase 
with 9^0 and then second order phase transition to the phase III. The phase transition 
between spiral and collinear phases is of the first and second orders at T < T* and T > T*, 
respectively, where T* is given by Eq. (I22l) . 
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IV. CONCLUSION 



We discuss in the present paper a nonfrustrated magnetoelectric ([2]) in magnetic field with 
spin-lattice coupling of the form ([T]) that is strong enough to produce spiral spin structure 
inducing ferroelectricity. Ground state energy is analyzed. We show that, in contrast to 
naive expectation, the onset of ferroelectricity takes place a.t H < Hd and Hd is lower 
than the saturation field ifc2 if constant (3 in Eq. ([1]) is not too large. The type of the phase 
transition between collinear paraelectric phase and spiral ferroelectric one depends on values 
of constants j given by Eqs. ([5]): the transition is of the second and first order if condition 
holds and does not hold, respectively. Moreover, if inequality (HM does not satisfy, 
there is only one critical field at which transition of the first order takes place from spiral 
phase to that in which all spins are parallel to the field. 

As a result of qualitative consideration using mean field approach we obtain phase dia- 
gram m H — T plane shown in Fig. [31 One goes on successively from plane (a) to (b) and (c) 
increasing (3. It should be noted also that the phase diagram remains qualitatively the same 
for a magnetoelectric with antiferromagnetic coupling among spins in xy-plane (cf. Eq. (I2D). 
In contrast, if the exchange coupling along z axis is ferromagnetic, the constant jz should be 
put equal to zero in the above consideration. As a consequence one leads to phase diagram 
shown in Fig. [3](c). 
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